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Abstract 

The symmetry properties under permutation of tomograms representing the 
states of a system of identical particles are studied. Starting from the action of 
the permutation group on the density matrix we define its action on the tomo- 
graphic probability distribution. Explicit calculations are performed in the case of 
the two-dimensional harmonic oscillator. 
PACS: 03.65.Wj, 05.30.-d, 02.20.-a 



1 Introduction 



In the conventional approach to quantum mechanics, that is the wave function descrip- 
tion ffl], the peculiarities of the quantum behaviour of a system of identical particles are 
encoded in the symmetry properties of the wave functions. The description in terms 
of density operators || also implies, for the density matrices associated to the states of 
identical particles, definite symmetry properties. Nevertheless both the descriptions are 
very different from the approach of classical statistical mechanics in terms of probability 
distributions defined on the phase space of the system. The wish to render the quantum 
picture closer to the classical one gave rise to hidden variables theories || and to the 
introduction of the Wigner quasi-distribution function on phase space J4| obeying to the 
Moyal equation ||. The same aim of describing quantum states in terms of a classical 
distribution of probability determined then the introduction of another class of quasi- 
distributions for the quantum states || [7|, |8|]. In 1969 Cahill and Glauber || considered a 
set of s-ordered quasi-distribution functions that further generalized the quasi-probability 
distribution functions previously introduced by Wigner, Husimi, Glauber and Sudarshan. 
They also showed that for each quantum state one can find some probability distribution 
which is determined by the density operator of the given state. But it was not clear 
whether or not such probability distribution determined uniquely the density operator, 
that is, whether or not the map was invertible. Only recently it was realized that quantum 
mechanics could be described completely in terms of this kind of probability distributions 
suitably defined for a random variable [Kj. In jTO] a consistent scheme has been proposed, 
the so-called probability representation, that results completely equivalent, in the sense 
of invertibility, to the ordinary formulation of quantum mechanics: the quantum states 
are described by a tomographic distribution of probability, or tomogram, also known as 
marginal distribution function (MDF), and the evolution of the system is described by an 
integro-differential equation for the MDF of generalized Fokker-Planck type. The prob- 
ability representation of quantum mechanics uses as a mathematical tool the symplectic 
tomography map [Tl[ of density operators onto quadrature probability distributions. For a 



general approach to tomograms and quasi-distributions of quantum states see also [T2|, [13 
and, with the inclusion of spin, |TJ|, [15], [16], |TT |. 

In the case of identical particles it is well known that they obey either Bose or Fermi 
statistics. As a result the wave function describing the system must be symmetrized or 
antisymmetrized respectively. Because the states of a physical system belong to a vector 
space and are, consequently, linearly superposable, it results quite easy to implement the 
two types of statistics in the ordinary quantum mechanics considering, for example, the 
wave functions of the system as a basis for a representation of the permutation group. 
On the other hand, in the probability representation of quantum states one needs to 
know what is the behaviour of the tomographic probability distributions for systems of 
identical particles. The aim of our article is the formulation of the symmetry properties 
of the tomograms for systems of identical particles. 
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In the following we use the relation between the MDF and the density matrix to 
construct the action of the permutation group on the marginal distribution functions. 
Because the set of the MDFs is not a vector space (see for example where a superpo- 



sition principle was formulated for both the sets of density matrices and tomograms) we 
obtain a realization, and not a representation, of the permutation group on the MDFs. 
In this way we can introduce the completely symmetrized and antisymmetrized marginal 
distribution function. 

The paper is organized as follows. In section 2 we review the realization of the permu- 
tation group on the set of density matrices and in section 3 we derive the corresponding 
realization on the set of MDFs. For the sake of clarity we work explicitly on the case 
of two identical particles and then we extend the results to n particles in section 4. In 
section 5 we describe a simple application and finally the concluding remarks. 



2 The permutation group of the density matrix 

Identical particles obey either Bose or Fermi statistics. This property implies two specific 
behaviours of the wave function of two identical particles: it must be symmetric in the case 
of Bose-particles, antisymmetric for Fermi-particles. These properties are more precisely 
described by means of the permutation group representation theory. To have a model let 
us concentrate on the system of two one dimensional particles whose positions are x± and 
x 2 respectively. The corresponding wave function, i[>(xi,x 2 ), can be decomposed into the 
sum of a symmetric and an antisymmetric function: 

1 1 

lf}{Xi,X2) = ~ (lp{Xl, X 2 ) + 1p{X 2 , Xi)) + -(lp{Xi,X2) ~ 1p{X2, Xi)) = 

i/j + (xi,x 2 ) + ij}-(xx,x 2 ) (2.1) 

This decomposition may be related to the irreducible representation of the permutation 
group G = (E,pi 2 ), where E is the identity of the group and p\ 2 represents the operation 
of permutation of the coordinates X\ and x 2 i.e. 

Ei/;(x 1 ,x 2 ) = ip(xi,x 2 ); 
Pi 2 ^{x u x 2 ) = ip(x 2l x 1 ). (2.2) 

The table of characters of the representations of the permutation group of two elements 
has the form 

E pl2 

1 1 (2.3) 
1 -1 

It means that the decomposition ( |2.1|) is connected with the irreducible representations 
Q2.2|) through the following formulas: 

Eip±(xi,x 2 ) = +lip + (x 1 ,x 2 ) 
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Pl2 r ip±(x 1 ,X 2 ) = ±llp+(Xi,X 2 ). 



(2.4) 



Thus the two functions ip + and realize a basis of the one-dimensional representation of 
the permutation group. We say that identical particles with symmetric wave function are 
described by means of the symmetric representation of G, while identical particles with 
antisymmetric wave function by means of the antisymmetric one. 

The density matrix p(x\, x 2 , x[, x' 2 ) of a pure state of two particles with wave function 
ip(xi,x 2 ) has the form 



p(xi,X2, x[, x' 2 ) = lj)(x\, X 2 )lf;*(x' 1 ,X 2 ) 



(2.5) 



thus we can extend the action of the permutation group to the set of density matrices 
because the coordinates X\,x 2 and be permuted independently. This means 

that the group G = G eg G which is the direct product of the permutation groups of two 
elements is the one related to the symmetry properties of the density matrix. G contains 
four elements: 



I = E®E, P 12 =p 12 ®E, P' 12 = E®p 12 

The action of the four elements on the density matrix 

Ip(x 1 ,x 2 ,x' 1 ,x 2 ) = 
P 12 p(x 1 ,x 2 ,x[,x 2 ) = 
P' 12 p(xx,x 2 ,x' 1 ,x' 2 ) = 

Using the multiplication table of the group G: 



II12 = Pl2 ®Pl2- 

is defined as follows: 



p(x 1 ,x 2 ,x' 1 ,x 2 ) 
p(x 2 , Xi, x 1: x 2 ) 
p(x 1 ,x 2 ,x 2 ,x[) 
p{x 2 , X\, x 2) X^J 



(2.6) 



(2.7) 
(2.8) 
(2.9) 
(2.10) 
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we may construct a realization through the action of G on the following distributions 

1 



f>2 



[ip + fip + P X2 p + P' 12 p 
l - [ip + fi P - P 12P - P' 12 p 



p- 



Pi = -(Ip-Up + P 12 p - P' 12 p 



Ip-Up-P 12 p + P' 12 p) . 



(2.11) 
(2.12) 
(2.13) 
(2.14) 



It is easy to see that such an action is given by: 
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Incidentally we note that p + and p_ are also obtained as: 

p + (x 1 ,x 2 ,x' 1 ,x' 2 ) = il> + (xi,x 2 )i[>l(x , 1 ,x' 2 ) = (2.15) 
- [(ipixx, x 2 )ip*(x[, x' 2 ) + ip(x 2 , 4) + x 2 )^*(x' 2 , x[) + ip(x 2 , xx)ijj*(x 2 , x[)] , 



p-(x 1 ,x 2 ,x[,x 2 ) =ip~(xi,x 2 )ip*_(x' 1 ,x 2 ) = (2.16) 
- [ip(xi, x 2 )ip*(x[, x 2 ) - i/)(x 2 , x 1 )^*(x[, x' 2 ) - ip(xi,x 2 )ip*(x 2 , x[) + ip(x 2 , x x )ijj*{x 2 , x[)] 

that is p + and p_ are the usual density matrices respectively associated to the usual 
symmetric and antisymmetric wave functions. (Here and in the whole paper everything 
is obtained for pure states but because of the linearity of the mixtures of pure states we 
can easily generalise to the density matrices of arbitrary mixed states of two particles.) 

Thus, what we have done so far is to extend the action of the permutation group 
(which we know for the space of the wave functions) to the space of density matrices 
(Eqs. ( |2 . 7| ) - ( pTTTJj ) ) by means of Eq. ( |2.5| ) which links in a simple way the density matrices 
to the wave functions. The natural generalization of this procedure, that is the definition 
of the action of the permutation group on the space of tomograms through the expression 
of tomograms in terms of density matrices, is not so simple since such a relation is an 
integral transform. This implies that the density matrices and the MDFs depend on 
different variables and, as we will see in the next section, it is not clear what is the action 
of the permutation group on the new variables. We follow then a different approach which 
we illustrate preliminarly for the space of the density matrices. What we are looking for 
is an expression of the symmetrized and antisymmetrized density matrices in integral 
form and, more precisely, we search for a realization of the permutation group in terms 
of integral kernels. Starting from Eqs. ( |2.7| ), ( |2.10| ) it is easy to realize that 



p± = ^{p(x 1 ,x 2 ,x' 1 ,x' 2 ) + p(xi,x 2 ,x' 2 ,x[) ± [p(x 2 ,x 1 ,x' 1 ,x' 2 ) + p(x 1 ,x 2 ,x' 2 ,x' 1 )}} 

= ^ J dyidy 2 dy[dy' 2 {5 yuXl 5 y2jX2 5 y > iiX > i 5 y > 2tX i 2 + 5 yuX2 5 y2 , Xl 5 y > i y 2 5 y2tX > i + 

± [Syi,x 2 Sy 2 ,x 1 S y ' 1>x ' 1 6y' 2iX > 2 + 5 yuXl 5 y2tX2 5 y > itX2 5 y2!X > i ] } ^(y a , j/ 2 , ^ , ^ ) , (2.17) 

with 5 yuXj = 8(yi- Xj), that is 
p±(xi, x 2 ; x[, x 2 ) = J K±(xi,x 2 ,x[,x 2 ;yi,y 2 ,y[,y' 2 )p(y 1 ,y 2 ;y[,y 2 )dy 1 dy 2 dy' 1 dy 2 , (2.18) 
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with the kernel 

K ± (x 1 ,x 2 , x[, x' 2 ; yi,y 2 , y[, y' 2 ) = ^ {<Wi ^2,^^,4 ^,4 + WAa.zAi,^^,^ 

,X2 fiy2,x 1 $y' 1 ,x'- i 3y' 2 ,x' 2 + Sy 1 ,x 1 ^y 2 ,X2^y' 1 ,x 2 ^y' 2 ,x' 1 '\ 19) 

Defining the permutation operator as 

^31 ijs—Jnfv*'!) "^2; • • • j X n ) f( y Xj 1 , Xj 2 , . . . , (2.20) 

we can write the kernel as: 

K±(X, x'\y, y') 2 = J^yi 1 ± P 12) ® (/ ± P^y^x^xjy'^y'^ (2.21) 

where (x, x'\y, y') 2 = (xi,x 2 ,x' 1 ,x' 2 ;yi,y 2 ,y' l ,y 2 ) and we assume that P, P' act on x, x' 
respectively. Analogously we may define the kernels Ki,K 2 associated to the density 
matrices pi,p 2 . In this way we obtain a realization of the permutation group in terms of 
integral kernels. This approach is absolutely equivalent to the one previously described 
but it is the one we need to implement the symmetry properties in the MDF framework. 



3 The permutation group of the MDF 



We are now in the position to extend the previous procedure to the tomographic proba- 
bility distribution, but we first need a brief review of the theory. 

The MDF of a random variable X is defined in f9H as the Fourier transform of the 
quantum characteristic function x(^) 



:< e ikx >; 



w{X, t) 



1 

27 



dke~ ikX < e lkX > 



(3.1) 



where X is the operator associated to X, and, for each observable O, < O >= Tr(pO) 
with p the time-dependent density operator. The MDF so defined is positive and normal- 
ized to unity, provided X is an observable 0. 

In Ref. |19| it is shown that by taking 



X = pq + up, 



(3.2) 



with q and p two conjugate variables and p, v real parameters labelling different reference 
frames in phase space, w(X, p,u,t) is normalized with respect to the X variable and 
there exists an invertible relation among the MDF and the density matrix. The variable 
X represents the position coordinate taking values in an ensemble of reference frames. 



Equation (3.1) may be rewritten in the more convenient form [20] 



w(X, p, z/, t) 



1ix\v\h 



p(Z',Z",t)exp 



,Z'-Z" 



x 



ii- 



z' + z" 



dZ'dZ" (3.3) 
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and its inverse is represented by 

p{X, X', t) = 7^J w{Y, l*,X- X') exp ^(y- t^^Y^j dn dY (3.4) 

It is important to note that, for (|3.3|) to be invertible, it is necessary that X be a coordinate 
variable taking values in an ensemble of phase spaces; in other words, the specific choice 
// = 1, v = or any other fixing of the parameters /i and v would not allow to reconstruct 
the density matrix. Hence, the MDF contains the same amount of information on a 
quantum state as the density matrix, only if Eq. fl3.2| ) is assumed. 

For Hamiltonians of the form 

an evolution equation governing the time dependence of the MDF is available (see ref. 
20fl for a simple derivation) 



d t w{X^,P,t) = ^d u + ~[V(-(d x )-%-^d x ) 

- V(-(d x )-% + ^d x )]^w(X,^i,,t) (3.6) 



where the inverse derivative is defined as 



{dxY 1 J f(Z)e^ x dZ = J f -^ e ^ x dZ. 



(3.7) 



The evolution equation, of generalized Fokker-Plank type, plays the role of the Schrodinger 
equation in the alternative scheme we are outlining. Its classical limit is easily seen to be 



^ ( » = { £i; + vV ' {- {m) 1 £) m] w(x ' t] 



(3.t 



where V is the derivative of the potential with respect to the argument. Equation (|3.8|) 
may be checked to be equivalent to Boltzmann equation for a classical distribution of 
probability f(q,p,t) , 

df p df dVdf , nn . 
at m oq oq op 
after performing the change of variables 

w{X,n,v,t) = — [ f(q,p,t)e lk{x - f " q - up) dk dq dp ; (3.10) 
2n J 

Hence, the classical and quantum evolution equations only differ by terms of higher order 
in ft. Moreover, for potentials quadratic in q, higher order terms cancel out and the 
quantum evolution equation coincides with the classical one. This leads to the remarkable 
result that there is no difference between the evolution of the distributions of probability 
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for quantum and classical observables, when the system is described by a Hamiltonian 
quadratic in positions and momenta. The generalization to N particles and eventually to 
field theory is straightforward and may be found in ref. [^0] . 

The scheme just outlined is selfconsistent and doesn't require at any step external 
structures such as the wave function or the density matrix. In this sense the probability 
description of quantum mechanics in terms of tomograms, once completely formulated, 
has to furnish a perfectly equivalent scheme to conventional ones |TU|. Thus, it appears 



quite natural to investigate the symmetry properties of the MDF of identical particles 
with respect to permutations. To this regard we discuss for simplicity the case of two 
particles, while a generalization will be exhibited in the next section. 

From eq. ( |3.3|) the MDF of two identical particles is represented by: 

w(i 1 ,p 1 ,u 1 ;i 2 ,p 2 ,u 2 ) = k 2 J D 2 {z;z')xe-^ ai+a ^p{z 1 ,z 2 ;z' 1 ,z' 2 ;t) (3.11) 
where we have defined 



ii) D n {z- z'-x-...- y ) = (nt 1 ^)(nt 1 ^)(n™=i^) • • • (n™ =1 %), 

m) k n = n™ =1 2^-y; 

and the time-dependence has been omitted. We note that because 

p{zu z 2] z[, 4) = z[)p 2 (z 2 , z 2 ) (3.12) 
the tomogram of a system of two particles factorizes as: 

6,A*2, y-i) = wi (£i, ^1,^1)^2(6^2,^2) (3.13) 

with 

Wi{ii,yLi,Pi) = * , J dzidz-e'^piizi,^. (3.14) 

Therefore it seems natural to identify the MDF relative to symmetrized and antisym- 
metrized states as 



w±(£,ii, v) 2 :=k 2 J D 2 {z- z'y-^^p^z') 



(3.15) 



with (£, /i, v) n = (£1, Pi, v\\ £ 2 , /i2, ^2; • • • in-, Pn-, v n ) (in the same way we may define wi and 
w 2 as the tomograms associated to p\ and p 2 . 

From (|2.18|) , which represents p± in terms of the integral kernel, we can write ( |3.15|) 
in the form 

w ± (i,p,u) 2 = k 2 J D^y-z^y^K^y'^z'^e-^+^piyly'h, (3.16) 



then, using the inverse formula ( |3.4|) for the case of two particles 

p(v\v'h = j^y 2 J D 2 (m, x)e^ + ^w(x, ™, y - y') 2 (3.17) 



with fa = Xi — mi Vl ~^ Vi we obtai 



am 



(2tt 

e -i[<*i+°>*}j\Pi+to\K ± (y, y; z, z') 2 w(x } m, n) 2 



= J D 2 (x,m,n)ij^^ J D 2 (y,y',z,z')5 niiyi ^ yl 5 n2>y2 ^ 2 

e- i[ai+a2] e i[ P 1+f52] K±(y, y'- z, z') 2 ^w(x, m, n) 2 . (3.18) 

Thus the tomographic probability distribution corresponding to a symmetric or antisym- 
metric state of two particles may be finally written as 



w±(£,jt/,i/) 2 = / D 2 (x, m, n)K± (£, /i, u\x, m, n) 2 w(x, m, n) 2 (3.19) 
with the kernel K± given by 



~ k f 

K±(£,n,u\x,m,ri)2 = J D 2 (y,y'; z, z')b n ^^ y >b n2 ^y 2 

x exj>[-i(a 1 + a 2 )+i(p l + (3 2 )}K ± {y,y'\z, z') 2 . (3.20) 



Also in this case it is straightforward but not particularly illuminating to derive the 
integral kernels Ki,K 2 , associated to w\ and w 2 , which complete the realization of the 
permutation group on the space of tomograms. 



4 Generalization to N particles 

The generalization to N particles is now straightforward: we define the kernel K±(y, y'; z, z') r 
as the following: 



K±(y,y';z, z') 



(4.1) 



where ep and epi are equal to 1 in the symmetric case, while representing, in the antysym- 
metric case, the sign of the permutations P and P' respectively. Therefore we obtain 



(2tt) t 



D n (y, y', z, z') IP ^ 



K±(y,y';z,,z') n (4.2) 



so that we find 



w±{£,iJ.,v) n = / D n (t,m,n)K±(£,[x,v\t,m,n) w(x,m,n) n . 



(4.3) 



We may introduce a sort of generalized Slater determinant to obtain K±, and, by 
means of fl4.2|), We define 



A±(y|x)„ = — 
n 



^Vi,xi 3yi,x 2 • • • 
fiy2,xi 3y2,X2 ■ • • 3y2,xn 

Vyn,xi ®y„,x2 ■ ■ ■ &y n ,x n 



(4.4) 



with the convention that for the symmetric case we take always the plus sign in computing 
the determinant. In this way we can write the kernel K±(y, y'\z, z') in a more conventional 
manner as: 

K±(y, y'; z, z') n : = A±(y\z) n A±(y'\z') n . (4.5) 

Summarizing, we have succeeded in realizing the symmetric and antisymmetric tomo- 
graphic probability distribution associated to a system of iV identical particles, eq. ( |4.3|) , 
in terms of the integral kernel ( |4.2| ) . This result fits into the selfconsistent scheme outlined 
at the beginning of this section in the sense of yielding a realization of the permutation 
group on the space of tomograms which doesn't require additional structures to be de- 
fined. In the forthcoming section we will see in concrete how it works on an example and 
we will check the invertibility retriwing the well known expressions for the symmetric and 
antysimmetric density matrices associated to a system of oscillators. 



5 Application 

Let us consider as an example the case of two independent harmonic oscillators. For this 



system the time-evolution equation ( |3.6| ) has been solved in ref. |2T[ and the solutions 
have the following expression: 

9 9 

e -vi-v£ 

W nm ( X ,fX,u) 2 = ^ rir2 \ n[m[2 n + rn H n(yi) H rn(y2) 

where Tj = e %t {\ij + ii/j), yj = Xj/\rj\, j = 1,2 and H n {y) and H m (y) are the Hermite 
polynomials of order n and m respectively. 

For the two-dimensional case the kernel may be obtained in a closed form by integrating 
Eq. ( p.20| ) over the z, z' and the y, y' variables. After some lengthy calculation we are 
able to write 

K±(£,fj-,vh = /(£i,A*i,fi; 6,^2,^2) + /(&, ^2,^2; 6, a*i, fi) 

± kK£i,A*i,^i; £2,^2,^2) +#(6,^2,^2; £1,^1,^1)] (5.2) 

with 
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/(&, 112, V2\ £l, 01, *l) = ¥ (^1 - ^ ) 5 - ) eX P (x X + X 2 - ^ - 



5 (6, ^2, ^21^1,^1,^1) 

u = {mi,ni,m 2 ,n 2 }; 



&2 ^1^2 
2tt 1^2^1-/^1^21 



^exp{iw T -A-w + J B-w + ii(xi + x 2 )}; 



/ l ^1^2 j M2 "1+^1^2 

» — /xa"l) 2(/i 2 ^l— pw) 

• /J2^1+/J1^2 j MlM2 
2(/i2^1— ^1^2) (Ml^2 -M2^l) 
















^ _/£2^1+/JO_^2_ 

' (jJblV2 -fJ.2Vl) 2(/i2Vl— ^1^2) 






^1^2 



\ 



M1M2 



B = { j ¥l&ZH2£l l M26-M16 ^ 1/1^2-^2^1 ^2gl-^l6 



2(^21/1-^1^2) (Mi^ -M21/1) 
L ^2^1-^16 

il\V1-ll2V\ ' " ^1 Vl-lllV\ ' ^2— M2"l ' Ml 

where t> T denotes the transposed vector and • is the usual matrix product, rows by 
columns. This is all what we need to obtain the symmetrized and antisymmetrized MDF 
corresponding to all the states of the two oscillators, although Eq. ( |3.19| ) is in general not 
easy to integrate. 

For the first non trivial case, i.e. the first excited state, u>±(£, /x, z/) 2 , we obtain, 
substituting Eq. Q and Eq. ([□]) for w 01 into Eq. |TT|) 



D 2 (x, m, n)K± (x, /x, u\x, m, n) 2 

6 



-1/1-1/2 



eg 



e ^i +l/ i ^1+ "2 



^v / (/"? + z/ i 2 )(/ u i + z/ 2 2 ) 

A*iA*2 + ^1^2 



v 7 /^ + ^ 2 



ir\rir 2 \n\m\2 n+m 
2 

2 



V 7 ^? + ^ 



± 2^6 



(/xf + ^^ + ^l) 



(5.3) 



As we can see, only the first two terms of the sum are probability distributions, whereas 
the last one is an interference term. 



For our edification we can check, using the inverse formula (|3.17|) , that we obtain the 
correct expression for the symmetrized and antisymmetrized density matrices associated 
to a system of two harmonic oscillators, that is 



p±(x,x') 2 = ^ exp | — ^(xl + x\ + x[ 2 + x' 2 2 ) \ [4x 2 x' 2 + Ax^ ± (4xix' 2 + 4x 2 x' 1 ) 



(5.4) 



6 Conclusions 



In summary, we have studied the permutational symmetry of tomograms of quantum 
states for systems of identical particles of bosonic and fermionic nature. We showed that 
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the tomographic probability distributions of identical particles in the probability repre- 
sentation of quantum mechanics are associated to the realization of the action of G®G, 
where G is the permutation group and £g> denotes the direct product. It is demonstrated 
that the tomographic probability distribution, that is the MDF is a sum of terms, some 
of them directly recognizable as tomograms, some others which are interference terms, 
that is not associated to a probability distribution by themselves, but only in the given 
combination. The main results of our analysis are summarized by Eqs. (4.2), (4.3), where 
tomograms of identical particles are explicitely derived in terms of an integral kernel. 
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